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^0 ' sessing a non-null Killing field It is known that for the vacuum case some of the basic field 
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equations are deducible from the others. It will be shown here how this result can be general- 
ized for the case of essentially arbitrary matter fields. The systematic study of the structure 
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^ ' of the fundamental field equations is carried out. In particular, the existence of geometrically 

, preferred reference systems is shown. Using local coordinates of this type two approaches are 
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presented resulting resolvent systems of partial differential equations for the basic field vari- 



^f-^ ' ables. Finally, the above results are applied for perfect fluid spacetimes describing possible 

, equilibrium configurations of relativistic dissipative fluids. 
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■ I. Introduction 



If one is given a complicated system of nonlinear partial differential equations to solve 
- as it happens frequently, for instance, in Einstein's gravitational theory - it is hard to see 
whether there exists any relationship between the equations or not. Sometimes the realization 
of certain type of connection might induce the introduction of an entirely new technique in 
solving the selected problem. This was the case, for example, when Cosgrove [1] gave a new 
formulation of fleld equations for stationary axisymmetric vacuum gravitational fields, or when, 
by a generalization of Cosgrove's approach, Fackerell and Kerr [2] derived a resolvent system 
of differential equations for the vacuum field equation of Einstein's theory for spacetimes with 
a single non-null Killing vector field. 
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In the first part of this paper we are going to show that the fundamental results the 
introduction of the new approach was based on in Refs. [1,2] can be generalized for spacetimes 
possessing a non-null Killing field with essentially arbitrary matter fields. Subsequently, the 
properties of the basic field equations are studied in the situation where the gradient of the 
norm of the Killing field and the twist of the Killing field are linearly independent. It is shown, 
for instance, that there exists a geometrically preferred vector field on the space of Killing orbits 
so that the basic field equations possess - in local coordinate systems adopted to this vector 
field - very simple form. In particular, a number of the relevant field variables and/or their 
partial derivatives with respect to the coordinate associated with the preferred vector field are 
found to be identically zero. By the application of the associated simplifications, two different 
approaches in deriving resolvent systems of partial differential equations for the basic field 
variables arc presented. The first is a general approach while the second one is a generalization, 
for particular matter fields, of the techniques applied for the study of stationary axisymmetric 
vacimm fields by Cosgrovc [1] . In the last part of this paper the application of both of these 
techniques for the case of perfect fluids possessing 4-velocity parallel to a timelike Killing field 
will be presented. 

II. The field equations 

In this section, first, we shall recall some of the notions and techniques of the formalism of 
general relativity developed for spacetimes possessing a non-null Killing vector field. Then, it 
will be shown that some of the field equations involved are always deducible from the others. 
Finally, as a direct application of this result the basic field equations will be reformulated - 
displaying the simplest form of the relevant equations corresponding to the possible subcases. 

Consider a smooth spacetimc, {M,gab), with a non-null Killing vector field, It is well 
known that for such a spacetime the formulation of the Einstein's theory can be simplified 
considerably by making use of a 3-dimensional formalism [3,4]. In particular, this is done as 
follows: Let S denote the space of Killing orbits of It is assumed here that <S can be given 
the structure of a 3-dimensional differentiable manifold so that the projection map, cf) : M ^ S, 
from M onto 5 is a smooth mapping [3]. This condition always holds locally, and for the 
case of a timelike Killing field in a chronological spacetime is shown to be satisfied globally [5] . 
Consider, now, the following three fields on M: the norm of the Killing field 



2 



(2.1) 



the twist of the KiUing field 

and the symmetric tensor field 



= Cabcd C^'V'^^'^, 



(2.2) 



hab = gab-V ^^a^b- (2.3) 

The images of these fields by the differential, 0*, of ^ give rise to tensor fields on the 3-space 
<S. For instance, ^^.hab is the natural induced metric on <S which is Lorentzian or Riemannian 
according to that the Killing field, is spacelike or timelike. (Hereafter we restrict our 
considerations to the 3-space «S so it should not cause a big confuse that the same notation will 
be used for the tensor fields living on S and for their natural 'pull backs' onto M.) 
Then the basic field equations are [4,5] 

Ral = lv-'D,DbV - ^V-\D,v)iDbV) + Iv-^LJaLOb - Kbiu^mCOn} + KKR'L, (2-4) 

Dla'^b] = -^abmnrKlt^'^^\ (2.5) 

D'^D^v = ^v-\Dmv){D^v) - v-^WmLo"" - 2R'ZrC, (2.6) 

D'^LUa = ^V-^UmD"'v, (2.7) 

(3) 

where i?^;, and Da denote the Ricci tensor and the covariant derivative operator associated 

(4) 

with hab, while, R^^ is supposed to be given in terms of the energy-momentum tensor. Tab, of 
the matter fields by virtue of Enstein's equations 

(4) 1 

Rab = 8n{Tab - -9abT). (2.8) 

Equations (2.4) - (2.7) relate the various typo of projections of the 4-Ricci tensor to tensor 
fields and their covariant derivatives living on the 3-space S. It is important that the entire 
geometrical content of Einstein's theory for a spacctime, (M, gab), with a non-null Killing vector 
field, can be uniquely represented by a 3-dimensional metric space, {S, hab), along with the 
fields V and uja satisfying the above set of field equations. Even more important that, to 
any 3-dimensional formulation, {(5, hab);v,uja}, of this type - up to gauge transformations - 
there exists a unique 4-dimensional spacetime, {M,gab), with a KiUing field, so that the 
projection map (j) : M ^ S reproduces the 3-dimensional formulation we started with. In fact, 
(2.7) is just the integrability condition ensuring that the 4-geometry can be recovered from the 
3-dimensional formulation [3,4]. 
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Note that no restrictions have been raised concerning the matter fields. In fact, what we re- 
ally need is that the matter fields be represented by tensor fields, ' 6 ^ 
and, a diffeomorphism invariant action be associated with them so that the energy-momentum 
tensor. Tab, and the Euler-Lagrange equations can be expressed in terms of appropriate varia- 
tions of this action. 

It is important to note that the invariance of Tab under the action of the isometry group 
associated with do not imply that the fields ^ ^i-,''^ ''"' are invariant. There are, for 
instance, exact solutions of the stationary vacuum Einstein-Maxwell field equations so that the 
electromagnetic fields are non-stationary [6] . On the other hand, whenever £5 ' ^ = 
for each value of i one might consider the unique decomposition of the fields ^(i,"^' 
into tensor fields which possess definite 'tangential' or 'perpendicular' character with regard to 
their free indices. These fields can be built up from tensorial products of and the pull 

backs of tensor fields ^ "i---"'" living on S. 

^(3) bi...b„ 

It is well known that equations (2.4) - (2.7) can be simplified by the introduction of the 
conformal metric hab defined as 

hab = evhab, (2.9) 

where e takes the value -|-1 (resp. —1) for spacelike (resp. timelike) Killing fields. Then (2.4) - 
(2.7) take the form 

Rab = lv-^{{DaV){DbV) + WaWb} + + SV-^habrC}R'Z, (2-10) 

D[a^b] = -ea6m„r<'"e, (2.11) 

D'^DaV = v-\{bmv){D^v) - u^^u^^} - 2ev-'RZrC, (2-12) 
D^LOa = 2v-'^oJmD"'v, (2.13) 
where Da and Rab are the covariant derivative operator and the Ricci tensor associated with 

hab- 

Although we are considering the set of basic field equations for spacetimes with a Killing 
vector field - in which case some simplification arise compared to the general case - the whole 
set of field equations is still rather complicated. For instance, equations (2.10) - (2.13) give 
rise - in local coordinates - to a system of coupled non-linear second order partial differential 
equations for the function v and the components of the tensor fields u)a,hab, ^ - . 

01... 

In fact, the situation is, in general, even worse because, in addition to (2.10) - (2.13), we have 
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to solve simultaneously the Euler - Lagrange equations which govern the evolution of matter 
fields in the spacetime. Note that these equations of motion, in general, couple to the above 
set of field equations increasing thereby the complexity of the whole problem. Therefore it is 
important to know what are the exact relationships between these equations. 

Now, we are going to show that in the formulation of Einstein's theory for spacetimes 

possessing a non-null Killing vector field the same type of simplification arises as for the case of 
stationary axisymmctric vacuum case realized by Cosgrove [1]. In particular, it can be shown 
that equations (2.10) and (2.11) are actually far more fundamental than (2.12) and (2.13). More 
precisely, by using (2.10) and (2.11) one can derive the following algebraic relationship 

[bkv) [b''ba,V-V-^{{bmV){b"'v) - W^W™} + 2ev-^d^^rC\ + [b^OJa - 2v-^iOmb"'v 

(2.14) 

Up to this point we have considered only the set of basic field equation (2.10) - (2.13) 
which are equivalent to Einstein's equations and an integrability condition. Remember that 

(4) 

the 4-Ricci tensor, i?^^ was assumed to be given in terms of the energy-momentum tensor. Tab- 
Thereby, the last term of the left hand side of (2.14), which is, in fact, ev~^hb"^V"'Tmn, cannot 
be put zero simply by referring to the 4-dimensional twice contracted Bianchi identity. It is 
known, however, that this term is identically zero whenever either the complete set of Einstein's 
equations or the Euler-Lagrange equations for the matter fields are satisfied. Since our aim is 
to derive a relationship between some of the relevant Einstein's equations to get rid of this 
term later it will be assumed that the equations of motion are satisfied by matter fields. [The 
author would like to say thank you to the unknown referee who pointed out the need for the 
clarification why the 4-dimensional Bianchi identity cannot be applied to set the third term of 
(2.14) to zero immediately.] 

The way one could get the relation (2.14) is the following: Substitute the right hand side 
of (2.10) for Rab into the following expression 

b-Rab - \bbR. (2.15) 
Then by using (2.11) a straightforward calculation yields that 

b'^Rab - ^bbR =l:v-4 {bbv) Ib'^bav - v-'{{bmv){b^v) - w^w-} + 2ev-'R'Zrc 
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+ COb 



b^OJa - 2v-^UJmb"'v 



(2.16) 



eV-\bbV){R'2^rn+V-^'^''D^a'^b] + [O'^pab - \bb{hr^Pmn) 



where Pab = {h'^K + ^v'^KbC^ORrnn- Since is a Killing field on M we get by (2.13) 

v-^uj'^D^a^b] = 2ev-V{VnU)R'p''e + ev-\bbv){R'^^rn- (2-17) 

(4) 

We also have, for instance, C^Rab = 0- Moreover, it can be shown by using the relationship 
between the covariant derivative operators Da and Da - with a tedious but straightforward 
calculation - that 

D'^Pab - iDbCh^^Pmn) = SV-^hb"" [V"Cn " ^ V„i?'''] - 2sV-^ hb"" nCm) Rp' e ■ (2-18) 

Now using (2. 16), (2. 17) and (2.18) we obtain 



WRab - \DbR =\^~A (Dbv) [o'^DaV - v-'{{Dmv){D^v) - w^w-} + 2ev-'R'ZrC 



(2.19) 



+ W6 



D^coa - 2v-'cOmD"'v\ + sy-' hb"' [V" R^^ - -VmR ] 



Since the tensor field Rab is just the Ricci tensor associated with the three metric, hab, - in 
virtue of the twice contracted Bianchi identity - we have that the left hand side of the previous 
equation is identically zero. This proves then that (2.14) holds identically. 

In the remaining part of this section we are going to study the consequences of the algebraic 
relation (2.14). We shall use the assumption that the Euler-Lagrange equations are satisfied for 
matter fields which implies in the case when they are derived from a diffeomorphism invariant 
action that 

V"T„6 = (2.20) 

so the third term of (2.14) is zero. Therefore, wo have that the relevant form of (2.14) says then 
that the above particular linear combination of the form fields, DaV and LUa, must vanish iden- 
tically. Correspondingly, there arc two subcases which have to be treated separately, namely, 
DaV and iVa might be either linearly independent or not. 

Whenever the two form fields and oJa are linearly independent only the trivial combi- 
nations of them can vanish identically. In this case (2.12) and (2.13) can be deduced from (2.10), 
(2.11) and the Euler-Lagrange equations. It is then sufficient to solve (2.10) and (2.11) along 
with the relevant equations of motion for matter fields since any solution of these equations will 
automatically satisfy (2.12) and (2.13) as well. 

Suppose now that the two form fields, DaV and uia, are linearly dependent. This might 
happen whenever one of them vanishes throughout or there exists a function, /, such that 

OJa = f- (DaV). (2.21) 
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a; Consider first the case of vanishing DaV, i.e., we suppose that v is constant throughout. 
Since we can introduce then a new KiUing field instead of by rescaling with an arbitrarily 
chosen constant factor we may assume here, without loss of generality, that v = e. Furthermore, 
for this case (2.14) implies that the relevant form of (2.13) is a consequence of (2.10), (2.11) 
and (2.20). Hence, the whole content of the basic field equations reduce to 

Rab = \w + + £/l„6rr}<l, (2-a-l) 

u;„a;" = -2eC„rr, (2.a.2) 

D{a^h\ = -ea6mnr<'"e- (2.^.3) 

/3; Suppose now that u>a = 0, i.e., is hypcrsurfacc orthogonal. Then (2.11) and (2.13) 
arc expected to hold, furthermore, the relevant form of (2.12) is simply a consequence of (2.14). 
Hence, the basic equations for the case under consideration simplify to [4] 

= \v~\b,v){b^v) + {K'^hl + et;-2/ia6rnCn- (2-/3) 

7; Finally, suppose that neither DaV nor iOa vanishes, and, there exists a function, /, such 
that (2.21) holds. Then the elimination of from (2.10) - (2.13) yields by using the above 
relationship 

Rah = lv-\l + f)ibaV){bt,v) + {h^K + eV-^habrC}RZ, (2.7-1) 

b'^b^v = v-\l - f){bmv){b^v) - 2ev-'RZrC, (2.7-2) 
{bj){b'^v) = + f){brr,v){b^v) + 2£<U"r], (2.7-3) 

{b^J){b,^v) = -CabmnrRp'^'e- (2.7-4) 

It can easily be checked that the relevant form of (2.14) implies that (2.7.2) is deducible from the 
Euler-Lagrange equations and (2.7.I), (2.7.3) and (2.7.4). Hence, for this last case, equations 
(2.21), (2.7.1), (2. 7. 3) and (2.7.4) display the entire content of the basic field equations. 

III. Geometrically preferred local coordinates 

In the remaining part of this paper we shall restrict our consideration to the case of in- 
dependent form fields, i.e., we suppose that {b^aV)oJb] 7^ on a subset <S of S. (The other 
possibility, when baV and iVb are linearly dependent, will be examined elsewhere.) According 
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to the results of the previous section to get a solution of the basic field equations, (2.10) - 
(2.13), it is sufficient to solve (2.10) and (2.11) along with the relevant set of Euler-Lagrange 
equations. In this section we are going to examine the properties of the fundamental equations 
(2.10) and (2.11). In particular, it will be shown that there exist geometrically preferred local 
coordinate systems in which these equations possess very simple form. 
We shall use the following shortened form of (2.10) and (2.11) 

Rab = \'V~'^{{DaV){DbV) + WaWb} + Pab, (3.1) 

-Cfa^b] = (Tab, (3.2) 

where 

Pab = {KK + ev-'habrC}RZ, (3-3) 



and 



(Tab = -tabmnC^Rp "C^- (3-4) 



Note that Pab is a symmetric while cfab an antisymmetric tensor field on S, both depending on 
the fields u, 

Since (£>[„«) w^j ^ on <S there exists a nowhere vanishing vector field. A;", there defined as 

r = e"-'"'{biv)oj^, (3.5) 

where Cabc denotes the 3-dimensional volume element associated with hab, i-e., Cabc = ^abcd^'^- 
Then the following hold 

CkV = 0, (3.6) 

k^Ua = 0, (3.7) 

k''{Rab-Pab) = 0, (3.8) 

and 

^k{Rab - Pab) = V~'^k''{aea>^b + (^ebOJa}- (3.9) 

Equations (3.6) - (3.8) are direct consequences of the definition of k". For (3.9) note that 

Ck{Rab - Pab) = k'De{Rab " Pab) + {Reb " Peb)Dak' + {Rae " Pae)Dbk^ ■ (3.10) 

However, according to (3.8) we have {R^b — peb)Dak^ = —k^Da{Reb — Peb), and so 

jCk{Rab - Pab) = k^{De{Rab " Pab) - Da{Reb - Peb) - Db{Rae - Pae)} ■ (3.11) 



Now, using (3,1), (3.6) and (3.7) we get 



which imply, along with (3.2) and the fact that Da is torsion free, that (3.9) holds. 
Note that whenever fc^cTeo is vanishing on <S, we have 

Ck{Rab-pab)=Q, (3.13) 
i.e., k°' is a collineation vector field of Rab—Pab- According to the definition of aab the contraction 

(4) 

k^aea is identically zero whenever there exist functions a, /? such that Rp "^^ = a^" + f3k'^. For 
the vacuum case, pab = (^ab = (or a = /3 = 0). Then (3.9) reduces to the well-known result 
that fc" is a Ricci collineation vector [7] . 

It is straightforward to check that (3.6) - (3.9) are satisfied not merely for k" but for any 
vector field possessing the form fk", where / is an arbitrary function on <S. Hereafter, the 
vector fields, k"" = fk"", which are defined with the use of a non- vanishing function, /, on <S will 
be referred as being geometrically preferred. 

Just like for the vacuum case (see Ref. [2]) one can introduce geometrically preferred local 
coordinate systems. Denote by any of the geometrically preferred vector fields and consider 
local coordinates, {x^ ,x'^ ,x^), adopted to /c", i.e., 

^" = (al^)"' = ^^-^'^ 

This type of coordinates can always be introduced (at least locally) on S. 

In such a local coordinate system, {x^,x'^,x^), equations (3.6) - (3.9) take the form 

W3 = 0, (3.16) 

^3/3 - P3/3 = 0, (3.17) 

and 

d 

Q^(^<^I3 ~ Pap) = ^"^{c^SaW/J + CT3f5U!a} (3.18) 

where P takes the values 1,2,3. Note that whenever one of the functions, p^/s (/3 = 1,2,3), 
does not vanish identically (3.17) gives algebraical relationship (s) between the variables v, hab, 



^ ai...afc g^jjj ^jjg derivatives of hab and ^ . Then we get that in such an 

'■'I bi...bk '■'> bi...bk 

adopted local coordinate system, {x^,x'^,x^), (3.1) is equivalent to (3.17) and 

Rab = ^v~'^{{dAv){dBv) +wawb} + PAB, (3.19) 

where Oav denotes the partial derivative of v with respect to the variable x"^, and the capital 
Latin indices take the values 1,2. 

It can be easily checked that in such a coordinate system (3.2) takes the form 

dlUJ2 — d^ijJl =(7l2 

(3.20) 

Observe that, whenever aab vanishes identically these equations imply that (at least locally) 
there exists a function, u) = oj{x^ jx"^), so that 

u)A = Oau). (3.21) 

Using these simplifications, in the next two sections two different methods in establishing 
resolvent systems of partial differential equations for the basic field variables will be presented. 

IV. General method 

This section is devoted to the introduction of a general approach to get a resolvent system 
of differential equations for the basic field variables. This approach is based on the following 
observation: It seems to be a general feature of the present formulation of Einstein's theory 
that (3.17) can be solved for the function v in many cases. Hereafter, we shall assume that 
the fields ^'f^,"^ "'"^^ are invariant under the action of the isometry group associated with 
thereby, we can use the fields ^ to represent the matter content, instead of 

them. Combining these two facts, hereafter we shall assume that the norm of the Killing field, 

V, can be given in terms of quantities derived from the induced 3-geometry, hab, and, possibly, 
from the tensor fields 'ipi^^"'^'"""^ ^ ^ , representing the matter fields. Correspondingly, we shall 
assume that there exists a function 

v = v{h^0, diKp, dsdX^; V'(,,""-"'"^^...^„, a^Vc,,"^-"-"^^...^^, ^J, (4.1) 

where the presence of second order partial derivatives of the fields ^ ^ indicates 

that the matter Lagrangian is supposed to contain at most second order partial derivatives of 
these fields and the Greek indices refer to components of tensor fields in geometrically preferred 
adopted local coordinates. 
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To start off note that (3.19) can be recast into the form 

Hab = v~'^{{dAv){dBv) +u;aojb}, 



(4.2) 



where 

Hab = 2{Rab - pab)- (4.3) 

It is important to emphasize that at each (cxphcit or impHcit) appearance of the function v in 
(4.2) the substitution of the right hand side of (4.1) is understood. Since we are dealing with 
the case of linearly independent form fields, i.e., {D[av)ijJb] ^ on <S, (4.2) can be shown to be 
equivalent to the following set of equations 

= HA.{d,v)-H,A{d,v) 

[det{HAB)]^ ^ 

and 

Hn{d2vf - 2Hi2idiv){d2v) + H22{divf - v'^det{HAB) = 0, (4.5) 

where the sign ambiguity of oja is indicated by the factor e (i.e., e = ±1) in (4.4). 

Using the definition, (3.3), of pab and (3.17) it can be checked easily that v depends on at 
most second order derivatives of the metric functions, /Iq/j, since only the terms Ras, enter (3.17). 
Therefore, with the assumption that at most second order covariant derivatives of the fields 
^ ai...o„ involved in the matter Lagrangian, we can conclude that (4.5) is at most a 

third order partial differential equation for the fields ha3 and V', i . Three additional 

partial differential equations restricting these fields have to be taken into consideration. These 
are derived by substituting the right hand side of (4.4) for uja into (3.20) and can be given as 
follows: 

{diH22){div) + H22{drd^v) - {diH,2){d2v) - 2H,2{d2div) - {d2H^2){div) 
+ {d2Hii){d2v) + H^iid2d2v) + ai (ln[detiHAB)]-^) (^^22(^1^;) - H^2{d2v)) (4.6) 
- d2(ln[det{HAB)]-'^) (H^2{div) - H^i{d2v)) = e[det{HAB)]^ ai2, 

{d3Hi2){div) - {d3Hii){d2v) + d3(ln[detiHABT^) [Hi2{div) - Hu{d2v)) 

= e[det{HAB)]^cr3i, 

{d3H22){div) - {d3H2i){d2v) + d3(ln[det{HAB)]-'^) (H22{div) - H^2{d2v)) 

= e[det{HAB)]^(T32- 

The first equation, (4.6), is a fourth order while the last two ones are third order non-linear 
partial differential equations. These equations along with (3.17), (4.5) and the relevant set of 
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Euler-Lagrange equations give rise to a resolvent system of field equations for the variables ha/3 
and t/Jc-i"^ a fl • Once one could get a solution of these field equations one can determine 
V via (4.1), moreover, wa can be given in virtue of (4.4). 

Clearly, the applicability of this approach strongly depends on the detailed functional form 
of V which was implicitly used throughout this section. For instance, the explicit form of 
the basic field equations, (3.17), (4.5)-(4,8), for the variables ha/3 and ^(j)"^' ""^ ^ can be 
examined only for particular matter fields separately. In section 7 we are going to give the 
functional form of v for perfect fiuid matter sources possessing 4-velocity parallel to a timelike 
Killing field and for particular equations of state. 

V. Generalization of Cosgrove's method 

In this section wc generalize the techniques developed originally for stationary axisymmctric 
vacuum fields for spacetimcs possessing a singe non-null Killing field with matter fields satisfying 
the additional conditions given below. More precisely, a slightly modified version of Cosgrove's 
approach will be established so as to derive from the basic set of field equations a resolvent 
system of differential equations for the basic variables. 

The two conditions are the following: 

Condition 5.1: The tensor field Uab vanishes throughout, i.e., ^[a-Rbje?^ = 0- 

Condition 5.2: The tensor field pab has the property that, in a geometrically preferred local 
coordinate system, its components, pab (^j B = 1,2), can be given exclusively 
in terms of the induced 3-geometry, hab- 

In particular. Condition 5.1 implies that (at least locally) there exists such a function w that 
= DaLO. Since we are dealing with the case of independent form fields - i.e., {D[av)iOb] 7^ ~ 
the functions v and w are then functionally independent. Condition 5.2 might be satisfied when 
(3.17) can be solved for v, moreover, (by using the relevant expression for v) one can eliminate 
thereby v and V'^)"^"'"'";, ^ from pab- Whenever both of the above conditions hold (3.19) 
can be recast into the form 

Hab = v-^{{dAv){dBv) + {dAUj){dBij)}, (5.1) 

where we used the expression (4.3) for Hab- Furthermore, due to Condition 5.2 the left hand 
side of (5.1) depends exclusively on the induced 3-metric while the right hand side of it depends 
merely on the functions v and w. Since v and uj are functions of a;^ and x^, equation (5.1) 
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shows that the same is true for the functions H^b even if, for instance, some of the components 
of hab may depend on x^. (Note that this property of the functions, Hab, is in fact a simple 
consequence of the general result (3.13).) Since v and iv are functionally independent we have 
that the functions Hab can be considered as the components of a non-singular Riemannian 
metric on a 2-dimensional manifold. Note that the right hand side of (5.1) is just the well-known 
representation of a Riemannian 2-metric in local coordinates {v, w) with Gaussian curvature —1. 
Hence, for the Gaussian curvature, , of the metric, Hab , 

= -1 (5.2) 

has to hold. This equation is, in fact, a fourth-order partial differential equation for the compo- 
nents of the tensor fields hat- For the case of linearly independent form fields under consideration 
(5.2) is the necessary and suHicient condition for the existence of functions v and iv satisfying 
(5.1). 

The outline of the proof of the above statement can be given as follows: Since we are 

considering the case of linearly independent form fields, (3.21) and (4.4) yield 

HA2idiv) - HiA{d2v) 
dAu; = ei 1 , (5.3) 

det{HAB)^ 

where the ambiguity in sign of uj is indicated by ei (i.e., ei = ±1). Substituting (5.3) into (5.1) 
with setting A, B = 2 and solving for div we obtain 

^ Hi2{d2v) + e2det{HAB)Hv^H22 - {d2vf]^ 
H22 

where H22 since otherwise {D[a'^)u)h] should vanish and €2 = ±1. Furthermore, the substi- 
tution of (5.4) into (5.3) yields 

d2Ui = eie2[v'^H22-{d2vf]i, (5.5) 

and 

det{HAB)Hd2v)-€2Hi2[v^H22-{d2vf]i 

diui = -ei ^ (5.6) 

-"22 

Equations (5.4) - (5.6) are equivalent to (5.1). The integrability condition, 828101 = 8182U), for 
the function cu can be shown [1] to give rise to the following Appel equation 

2H22{d2d2U) - AH22{d2Uf - (52i?22)(52«7) + HI2 + $[i?22 - 4(52f/)2]* = 0, (5.7) 

where U = ^ln{sv) and 

$ = ^£2 • det{HAB)-^{-2H22id2H2l)+H2lid2H22)+H22{dlH22)}. (5.8) 
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Utilizing Cosgrove's substitution (see Ref. [1,2]) 

d2U =-{H22)^M{l + M^)-\ 



(5.9) 



we obtain from (5.4) and (5.7) the following pair of Riccati equations 



OaM = Xa + 2YaM + ZaM'^. 



(5.10) 



Here the functions Xa,Ya and Za are defined as 



Xa 



Met{HAB)^H22 





(5.12) 



and 



Za = Xa - Ha2{H22) ^ ■ 



(5.13) 



The integrability conditions for the simultaneous set of Riccati equations, (5.10), reduce to a 
single condition [1,2], which, not unexpectedly, may be put into the form of (5.2). 

Summarizing the results of this section we can say the following: To get a resolution 



simultaneous Riccati equations, (5.10), can be used to determine the function v via (5.4) and 
(5.9). Afterwards, (5.3) can be applied to construct the function u. A detailed discussion 
about the resolution of the corresponding problems for the vacuum case - particularly, about 
the solutions of Riccati equations of the above type - can be found in Ref. [1]. Finally, 



components of tensor fields representing the matter content. 

It is worth mentioning that (5.1) inherits a remarkable feature of the corresponding equa- 
tion given for the vacuum case noticed by Gcroch [3]. Namely, this equation is invariant under 
the action of an 5*^(2,11) transformation. Two of the relevant parameters are associated with 
gauge transformations but there exists a one-parameter subclass of 'effective' 5*^(2,11) trans- 
formations yielding new solutions from known ones. In particular, by starting with a particular 
solution, {vo,u)o), of (5.1) associated with a fixed set of functions Hab one can generate a 
one-parameter family of solutions, {vT,0Jr) to this equation. More precisely, one can show by 
a straightforward modification of the proof of Theorem 1. of Ref. [1] that for fixed functions 
Hab satisfying (5.2) the full set of solutions of (5.1) - apart from those related to gauge trans- 
formations of the spacetime, {M,gab), - is generated from the particular solution, {vo,cuo), by 



of the basic field variables we have to solve first (5.2) for haj3- 



Then the solutions of the 



using these functions - w, w and ha^ - the Euler-Lagrange equations have to be solved for the 
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the transformation 



(5.14) 



= -, : r^— g • o • (5.15) 



{cosT — LOosinry + vfysin^r 

{siriT + (jJocost){cost — ajosinr) — v^sinrcosT 
{cost — iVosinr)'^ + VqSivP't 



There is, however, a significant difference between the vacuum case and the case under 
consideration. Namely, for the case of vacuum the relevant form of (5.1) is the only field equation 
to be solved while for the general case with matter the basic field variables have to satisfy, 
beside (5.1), both (3.17) and the relevant set of Euler-Lagrange equations, as well. Therefore, 
one would expect that there is no matter field so that the above transformation can be applied. 
Nevertheless, there exists such a matter field (see section 7) where certain restrictions on the 
basic field variables (associated with the matter content) can ensure the applicability of the 
transformation (5.14)-(5.15), and, consequently, one may generate new solutions of Einstein's 
equations. In particular, this transformation was used to derive a number of new perfect fluid 
solutions from known ones [8,9]. 

VI. Perfect fluids 

In this section some of the basic notions and results in connection with perfect fluids will 
be recalled and some of the consequences of the presence of Killing fields in the spacetimes will 
be discussed. 

Consider a perfect fluid with mass density, p, and pressure, P, (both quantities measured 
in the rest frame of the fluid), furthermore, with 4- velocity u", where u°'Ua = — 1. (Note that 
the tensor fields ^ ^i)""^'"""^ ^ on M for the present case are the fields p, P and u".) The 

energy-momentum tensor is given as 

Tab = PUaUb + P{gab + UaUb), (6.1) 

furthermore, the Euler - Lagrange equations are 

«"Va/9+(/9 + P)V"«a = 0, (6.2) 
[p + P)w"VaW6 + [gab + UaUb)^ P = 0. (6.3) 

It is known that for perfect fluid sources these equations are equivalent to the 'integrability' 
condition of Einstein's equation 

V"T„6 = 0. (6.4) 
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In particular, (6.2) and (6.3) are equivalent to the 'parallel to u"' and the 'orthogonal to u"' 
projections of (6.4), respectively. Thereby, it is usual in the formulation of Einstein's theory for 
spacetimes with perfect fluids to postulate merely the form of the energy- momentum tensor. Tab, 
and solve Einstein's equations since the equations of motion for the fluid then are automatically 
satisfied. We have chosen, however, a somewhat reversed approach here. In section 2 it was 
assumed that Euler-Lagrange equations are satisfied (which implies for the present case that 
VTa(, = 0) and this condition was used to show that some of the basic field equations are 
deducible from the others. It is important to emphasize that we earn more than we loose 
by replacing the two basic field equations, (2.11) and (2.12), by Euler-Lagrange equations. 
Equations (2.11) and (2.12) are second order partial diff'erential equations while the above 
Euler-Lagrange equations are first order ones for perfect fluid. 

Consider now the consequences of the presence of a Killing field, £,'^, for perfect fluid matter 
sources. First of all, 

C^Tab = 0. (6.5) 

Again, by the presence of a preferred vector field, u", one might consider the unique decom- 
position of C^Tab into symmetric tensor fields so that each of these tensor fields has definite 
'tangential' or 'perpendicular' character with regard to their free indices. Since C(Tab vanishes 
all of these projections must vanish, as well. Thereby {C(Tab)u'^u'' = which gives that 

C^p = 0. (6.6) 

Then (£^T„6)M"7r^ = yields that 

£{u" = 0, (6.7) 

or p + P = 0. (Note, however, when the equation of state is chosen to be p -I- P = then 

the energy-momentum tensor is of the form Tat = PQab, and (6.4) implies that P is constant 
throughout. This is precisely the case of vaeimm fields with non-zero cosmological constant 
so it seems to be reasonable to assume that p -|- P is not identically zero, and hereafter we do 
that.) Finally, from {C^Tab)'!T°' e'^'' f = we get 

C^P = (6.8) 

throughout, where the projector, Tr'^b, is defined to be 7r°'b = S°'b — u°'Ub. All in all, each of 
the physical quantities related to the perfect fluid are invariant under the action of the isom- 
etry group associated with Consequently, for a general perfect fluid spacetime possessing 
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a non-null Killing field we can use, without loss of generality, instead of the fields p,P, 

given on M the fields p,P, w,, = w^^a, u\ = h\u^ {''^^u)""""'" b,...C^ 

defined on S. 

Determine now the relevant form of Pab and Gab- According to (2.8) and (6.1) we have 

Ral = [{P + P)UaUb+\{p- P)gab] , (6.9) 

furthermore, by the definition of pab and (Tab 



Pab = 8n 

and 



(p + P){{ha"'Um)ihb^Un) + SV-\u%fhab} + eV-\p - P)hab , (6.10) 



(Tab = -87re„6,n„r«"(p + ^)(w^Ce), (6.11) 



hold. 

For simplicity, one may restrict ones considerations to the case of vanishing (Tab- Equation 
(6.11) implies that aab = whenever either of the following hold: Ug^" = 0, or ^["^u**! = 
(or p + P = but this case has been excluded earlier). Thereby, we can say that aab = 
throughout if and only if either the 4- velocity of the fluid, u", is parallel to the Killing field, 
which means that the spacetime is stationary and 

w'' = (-v)-i^\ (6.12) 

or 

M"^a = 0, (6.13) 

which might be the case whenever the Killing field, is spacelike. For both of these cases 
Condition 5.1 holds which implies that there exists (at least locally) a function oj such that 

tOa = Dato. (6.14) 

Let us consider the following particular case of perfect fluid sources: There arc two com- 
muting Killing fields, i;"^^ (A = 1, 2), on the spacetime and the 4-velocity of the fluid, u"', can 
be given as a linear combination of these Killing fields 

«« = >l(^«,+B^«J. (6.15) 

Then with linearly independent Killing fields (6.7) and (6.15) yield that the functions A and B 
satisfy 

Ci^^A = C^^^^B = (A = 1,2). (6.16) 
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Now, applying (6.7) and (6.8) for the Killing fields, , and using (6.15) we get that 

LuP = CuP = 0, (6.17) 

and, the equations of motion, (6.2) and (6.3), reduce to 

(p + P)V<'«„ = 0, (6.18) 

{p + P)u''VaUb + VbP = Q. (6.19) 

One extracts from (6.15) - (6.16) that the fluid is expansion free, i.e., VaW" = throughout. 
Thereby, (6.18) holds identically. Furthermore, a straightforward calculation yields that 

«-VaW6 = -\A^[Vb{-A-^) - '-' Vfeg}. (6.20) 

which along with (6.19) and (6.20) gives that 



1 



2^ 



VaP + ^ip + P) V„(ln^-2) - ' VaB = 0. (6.21) 



dB 



As it was argued in Ref. [10], (6.21) implies that P = P{A, B) and p = p{A^ B) even if A and 
B are functionally dependent or constant. Furthermore, since the 4- velocity - given by (6.15) 
- is a unit timelike vector we have 

A-' = -{(Cr„Ca,a) + 26(Cr„^(.,a) + (6-22) 

which along with (6.21) (and the above conclusion) gives that the equation of state must be of 
the form 

P = P{P)- (6.23) 

The remained Euler - Lagrange equation, (6.21), simplifies further whenever ^^'^"^ — ^VoB = 0, 
i.e., 

V„B = or ^ii^^=0. (6.24) 

The case VqB = is that of a 'rigid fluid', i.e., the 4-vclocity, is parallel to the timelike 
Killing field = + B^'^ . It is important to emphasize that equations (6.15) - (6.23) along 
with their consequences hold (with B = 0) without any alteration even if the spacetime admits 
only a single timelike Killing field, ^" = parallel to the 4-vclocity of the fluid, m". 

The other possibility, ^^^gg''^ = 0, along with (6.22) gives that B = - i.e., the 

«(2)«(2)/ 

4- velocity of the fluid, u", is orthogonal to which, therefore, must be a spacelike Killing 

(2) 

field. 
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Note that for both of these cases not merely the Euler-Lagrange equations are simplified 
but, in accordance with this fact, the potential space associated with the Lagrangian of this 
particular case of 'gravity plus perfect fluid' system admits a symmetry [10]. 

Moreover, equations (6.21), (6.23) and (6.24) yield then that 

rP 



(6.25) 



where Aq and Pq s^re constants of the integration. Consequently, whenever the 4- velocity of 
the fluid, m", is either parallel to a Killing fleld, or spanned by two commuting Killing fields, 
as in (6.15), with B = — y'"' , and, the equation of state, p = p{P), is known then the 

«(2)^(2)/ 

function A~'^ = A~'^{P) or P = P{A~'^) can be determined via (6.25). Note that the function 
A~^ possesses the form 

if mI^^^I = 0- 

t ac n (6-26) 



where 

= -KA^^eMiA^.f) + ii'^^A^^Hf. (6.27) 

Note that the function W has the following simple geometrical meaning. In canonical Weyl 
coordinates, (p, z, (j)), the 3-metric, hab, can be given as 

hocfi = diag{Exp{2j), Exp{2j), -W^}, (6.28) 

where 7 and W are functions of the coordinates (p, z) [4]. 

VII. Perfect fluids with 4-velocity parallel to a Killing field 

In this section we shall apply the results of the previous sections for perfect fluid spacctimes 
possessing a timclikc Killing field, parallel to the 4- velocity of the fluid, Such a fluid 
has expansion- and shear-free flow, i.e., it is 'rigid'. Thereby one might ask whether there 
exists any physically realistic situation in which such a model can be applied. However, it was 
shown by Geroch and Lindblom [11] that in a generic theory of relativistic dissipative fluids 
the equilibrium states are perfect fluid states. Furthermore, they showed that for these perfect 
fluids - which represent the equilibrium configurations of dissipative relativistic fluids - the 
4-velocity is parallel to a Killing field [11]. Therefore, the model we are dealing with in this 
section has to have physical relevance, and, in fact, it is the needed one as long as we are looking 
for a faithful description of possible equilibrium configurations of relativistic dissipative fluids. 
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First the applicability of the generalization of Cosgrove's method then the general approach 
will be considered. Clearly, for this type of perfect fluids, Condition 5. 1 is satisfied and we show 
that Condition 5.2 holds, as well. Now, since ul"^*"! = 0, the Killing field, is timelike so e 
takes the value —1. Furthermore, (6.6) yields that 

Pab = lewv-^Phab. (7.1) 

The relevant form of (3.17) 

= 16TTV-'^Ph3f3, (7.2) 

can then be solved for v. Since we have a non-vanishing spacelike vector field, fc", and the 3- 
metric , hai3, is non-singular, /133 cannot vanish. Whenever there is another non- vanishing one 
among the functions, /13/3, then (7.2) gives rise to an algebraical restriction on the components 
R30 of the Ricci tensor associated with hat- We obtain from (7.1) and (7.2) 

PAB = ^^hAB, (7.3) 
^33 

which means that Condition 5.2 is satisfied. Furthermore, this equation yields, along with (4.3), 

Hab = 2{Rab - ^Hab)- (7.4) 
/133 

For the particular case under consideration the functions Hab depend exclusively on the 
induced 3-metric, hab, and the relevant form of (5.2) is, in fact, a fourth order partial diff'erential 
equation for the components of hab- It is striking to what an extent the corresponding basic 
field equations are similar in structure to the vacuum counterparts. Turning back to the main 
issue, note that the functions v and u) can be determined by virtue of (5.3)-(5.13). Finally, 
the pressure, P, can be determined by (7.2) and the mass density, p, by the Euler-Lagrange 
equations, (6.21). 

After solving (5.2) and fixing the functions Hab we may ask for the conditions under 
which the transformation (5.14)-(5.15) yields new solutions of the basic field equations. The 
two equations to be solved are, for the present case, (7.2) and (6.21). It is straightforward to 
check that by choosing Pr to satisiy the equation PtV~^ = PqVq^ and deriving Pt - for each 
pair of the functions Vr and Pr - in virtue of (6.21) we get a one-parameter family of solutions of 
the basic field equations. Note that the invariance properties of the Lagrangian of electrically 
charged rigid perfect fluids was studied earlier by Kramer, Neugebauer and Stephani [4,12]. 
However, their considerations were restricted to the case of a static spacetime exclusively. 
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Thereby the transformation transparented by eqs. (5.14) - (5.15) is new for perfect fluids 
although it is a straightforward generaHzation of a known algorithm. Further analyzes related 
to this transformation and derivation of new perfect solutions by making use of it can be found 
inRef. [8,9]. 

An interesting subcase of these perfect fluid spacetimes, discovered many years ago by 
Ehlers [12,4], is the case of vanishing pressure, i.e., of a stationary spacetime with dust possessing 
4-velocity, u", parallel to a timelike Killing fleld, According to (3.17), (3.18) and (7.1) the 
field equations are then the same as for the vacuum case. Note, however, that (6.21) implies 
then V = const. The value of v may be chosen throughout to be -1, so the fleld equations 
simplify to 

Rs(3 = 0, (7.5) 

RaB = ^{dAUj){dBUj). (7.6) 

For the energy density, p, we obtain from (2. a. 2), (6.9) and (6.12) the constraint 

p = ^id^u;)idAu;). (7.7) 

Accordingly, a stationary dust {sd) solution - represented by v''"''^ = —1, w'^**' = ln[w'°'''] and 
^ab^ = /i^J ' - can be assigned to every static vacuum (sv) solution - given in terms of w''*^' and 
^ab ~ where the energy density of the dust satisfies the constraint (7.7) [13,4]. 

Although the above general method can be applied in a straightforward way to get solutions 
of Einstein's equations for the selected perfect fiuid source there is an unfavorable aspect of 
this method. Namely, the most significant physical quantities - the mass density, p, and the 
pressure, P, characterizing the possible physical states of the fluid - can be determined only at 
the very end of the entire process in terms of the function v and the 3-geometry, hat- Therefore 
the equation of state, p = p{P), has to bo in accordance of the corresponding constraints, which 
implies that it cannot be chosen freely. The general approach, introduced in section 4, can be 
used to cure this problem but the price we have to pay is the appearance of extra non-linearities. 

The general method was developed to ensure more control on the physical properties of 
matter flelds in solving the relevant set of field equations. The significance of the seemingly 
technical differences between the two methods can be transparented for the examined perfect 
fiuid sources as follows: Rcnicimber that the basic point in the general approach was the speci- 
fication of the functional form of v in terms of the 3-geometry and tensor fields representing the 
matter content. Moreover, for the present situation the basic field equations are (5.1) - where 
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Hab is given by (7.4) -, (7.2) and (6.25). Note that (6.25) - which is, indeed, the integrated 
form of the Euler-Lagrange equation - gives now the functional relation between the functions 
V, p and P as 



r dP' 1 

.^P)=.,.e.,[-2j^-^^]. (7.8) 

With the help of (7.2) and (7.8) it can be shown that the function v can be given in terms of 
the 3-geometry exclusively. Note that (7.8) is independent of the 3-geometry and, indeed, it is 
the only equation where one can control the physical properties of the solution describing the 
fluid by the substitution of a physically realistic equation of state. Consider, for instance, the 
case of polytrope equation of state, i.e., 

p = cP«, (7.9) 
where a G (0, 1) U (1, oo) and c > 0. Using then (7.8) and (7.9) one gets 

- pa-l I -| 

V = vo . (7.10) 



After solving (7.2) for v and substituting the resulted function into (5.1) one might attempt 
to get solutions of the yielded equation, where now Hab is given by (7.4). Remember that 
for the present case of perfect fluids aab is identically zero thereby one has to solve merely 
the relevant form of (4.5) and (4.6) along with the possible two equations involved by (7.2) 
[Note that the last two equations might give further algebraical relationships between certain 
components of the 3-metric and the 3-Ricci tensor.] Although the derivation of the equations 
is straightforward the appearance of extra non-linearities - related to the functional form of v 
- are frightening. Note, however, that whenever one is able to find solutions of these equations 
the physical relevance of the solutions is automatically assured. 

VIII. Final remarks 

A new formulation of Einstein's equations for spacetimes admitting a non-null Killing vec- 
tor field and arbitrary matter field was given in this paper. First it was shown that some of the 
basic field equations are always deducible from the others. Then the existence of a geometri- 
cally preferred vector field and related coordinate systems were shown. Based on the associated 
simplifications, two methods were presented obtaining systems of partial differential equations 
for the basic variables associated with the spacetime geometry and with the matter content. 
Both of the developed approaches were applied for perfect fluid spacetimes which describe equi- 
librium configurations of relativistic dissipative fluids. The symmetry properties of the relevant 
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equations and differences of the two approaches were studied. It was shown, furthermore, that 
the techniques which were developed by Cosgrove [1] for the vacuum stationary axisymmetric 
problem can be generalized straightforwardly for these perfect fluid spacetimes despite the fact 
that in our examinations we assumed merely the existence of a single timelike Killing field. 

It is worth emphasizing that the general results of this paper, given in details in sections 
2-5, are valid for any spacetime in Einstein's theory which possesses a non-null Killing field 
and essentially arbitrary matter fields. Thereby, it would deserve further studies to find out 
how to apply these results for even more interesting situations in which time dependence may 
occur and/or different types of matter fields are present. 
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